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Abstract. The search for muon to electron conversion with a sensitivity of the order 
10^^^ requires a several order of magnitude increase in muon intensity and a high 
resolution, cr ~ 0.1 MeV/c of the electron's momentum. 

We present results of a pattern recognition and track momentum reconstruction 
algorithm that relies on a Kalman filter approach. Background from captured protons, 
neutrons, photons and from muon decay in orbit were generated by GEANT. The 
effective average straw tube background rate was 800 kHz. 

The pattern recognition proceeds in two stages. In the first, simple considerations 
using only straw tube center coordinates, without drift time information, were applied 
to reduce the background to a manageable level. Then the drift time information is 
incorporated and a Deterministic Annealing Filter applied to reach the final level of 
background suppression and to provide a starting point for the track momentum recon- 
struction using the Kalman filter. This procedure reduces the simulated background 
by a factor 800 with small, (2. 7%), losses in real tracker hits. 

The momentum resolution of the tracker is cr = 0.12 MeV/c and the acceptance for 
muon conversion events with momentum above 103.6 MeV/c is 22%. These numbers 
do not differ significantly from the values obtained without background. 

The expected number of events from muon decay in orbit (main background) in 
which the decay electron has momentum greater than 100 MeV/c is 0.3, compared to 
6.5 /ie - conversion events above the same threshold for i?^e 10~^^. 



I INTRODUCTION 

The observation of fie - conversion would provide the first direct evidence for 
lepton flavor violation in charged lepton sector and require new physics, beyond 
the Standard Model (see [1] and references therein). Lepton flavor is not conserved 
in neutrino oscillations but the modifications to the Standard Model to include the 
small neutrino masses do not lead to an appreciable rate for /ze - conversion 

fi- + N e~ + N . 

This process violates the lepton flavor numbers, Lg and L^, but conserves the 
total lepton number. The signature of the process is very clear: a single monochro- 
matic electron in the final state with the energy close to the muon mass: 

Eg "f^ii -B^t -^rec 



^) Permanent address: Institute for Nuclear Research, 60-th Oct. pr. 7a, Moscow 117312, Russia 



where is a muon mass, -B^ is a binding energy of the Is muonic atom, Erec is a 
nuclear recoil energy. 

The SINDRUM II collaboration at PSI has carried out a program of experiments 
to search for fie - conversion in various nuclei. They find that at 90% CL the upper 
limit for the reaction + Ti — > e~ + Ti is 6.1 x 10~^^ [2]. In this experiment 
muons were accumulated at a rate of 10^ ii~ /sec. According to preliminary results 
[3] for the reaction /j," + Au — > e~ + Au a single event sensitivity is 3.3 x 10~^'^ 

In [4] an idea of increasing of muon beam intensity by a few orders of magnitude 
up to 10^^ ii~ / sec based on the solenoid-capture scheme was discussed and MELC 
experiment was proposed [5] with a goal to reach a sensitivity of the order 10"-^^. 

A new jie - conversion experiment MECO (Muon Electron Conversion) E-940 
[6], exploiting the idea of the solenoid-capture scheme, is under preparation at 
BNL. MECO aims to search for /i~ + A/ — > + Al with a single event sensitivity 
2x10^^^. It will use a new high-intensity pulsed muon beam, which could yield 
about 10^^//~/sec stopped in a target. 

Also the PRIME (PRISM MuE conversion) working group at KEK expressed an 
interest [7] to carry out a search for lepton flavor violation in /ie - conversion in a 
muonic atom at a sensitivity of 10~^^ using a proposed high intensity pure muon 
source of 10^^ — . 

In this article we develop a pattern recognition and track reconstruction proce- 
dure based on the Kalman filter technique for a transverse version of tracker for 
MECO experiment. 

This paper is organized as follows: In section 2 the transverse tracker is described 
and the advantages of the tracker discussed. In Section 3 possible backgrounds are 
discussed and the procedure of background simulation briefiy explained. In Section 
4 the pattern recognition procedure is developed. At this stage a deterministic an- 
nealing filter (DAF) is applied to make a final background suppression and provide 
a starting point for track momentum reconstruction by the Kalman filter tech- 
nique. Section 5 describes the procedure of momentum reconstruction based on 
the Kalman filter. Results of the pattern recognition and momentum reconstruc- 
tion are presented in this section. A brief summary outlook can be found in Section 
6. The appendices provide a more detailed look at Kalman filter and deterministic 
annealing filter. Also tracker resolution is discussed in Appendix. 

II THE TRANSVERSE TRACKER DESCRIPTION 

The goal of the MECO tracker is to detect the electron from /ie - conversion with 
large acceptance and measure its momentum with high resolution (cr ~ 0.1 MeV/c). 
The tracker is located in a uniform IT magnetic field. The minimal tracker length 
is defined by the following requirements: 

Background electrons with an energy around 105 MeV are produced by cosmic 
rays in the wall between the transport and detector solenoids. These electrons 
cannot have a pitch angle 6 in the tracker greater than 45°, due to the adiabatic 
character of charge particle movement. To suppress the cosmic ray background we 
require that measured pitch angle 9min for signal events to be more than 45°. 

To get more redundancy for signal events we require that the measured trajectory 
should have two full turns. This requirement sets limits on the minimal tracker 
length, which is expressed as: 

^Tracfeer = 4 • TT • 1/2.998 • P ■ cos{9min) / B ~ 4.19 • 105 • 0.707/1.0 ~ 310 cm 

o 



where P is a muon conversion momentum (MeV/c) and B is a magnetic field (1 
Tesla) in the tracker region. 

The tracker consists of 18 modules spaced 17 cm apart. A module consists of 6 
planes, each turned at 30° and shifted 2.5 cm relative to the previous one (Figure 
1). A plane consists of two trapezoidal chambers of width 30 cm and lengths 70 
to 130 cm in an up and down configuration. The chamber coordinate systems 
are defined by a rotation angle proportional to 30° giving an effective "stereo" of 
crossed directions for 12 different views. 

A protective sheath is used to suppress background Compton electrons from a 
chamber frame. The protective sheath makes nonsensitive a small region of the 
anode wire near the frame. 

The chamber consists of one layer of straw tubes (60 straws) of 5 mm diameter, 
and length varying from 70 cm to 130 cm. The total number of chambers is 216. 
The chamber sensitive area starts from 38 cm radius. The straws are assumed to 
have wall thickness 15 //m and are constructed of kapton. The total thickness of 
each chamber is 9 mg/cm^. The total number of the tracker straws is 12960. The 
tracker length is 302 cm. A signal from the straw anode wire will be used to get 
drift time. 

The Al target is tapered in the downstream direction, with 5 cm disk spacing 
and radii from 8.3 cm to 6.53 cm. The target is placed in the graded portion of 
the DS magnetic field, with the first disk at 1.75 T and the last at 1.3 T. Protons 
from the muon capture are absorbed in three concentric polyethylene absorbers: a 
conical tube of dimensions Ri = 46 cm, R2 = 70 cm, L = 260 cm and a tube (R ~ 
70 cm, L = 200 cm) both have thickness 3 mm, and a tube of smaller radius (R = 
36 cm, L = 235 cm) of thickness 0.5 mm which are placed before the tracker. 




The transverse tracker is to be compared with a longitudinal version that MECO 
is also considering. The longitudinal tracker [6] consists of an octagonal array of 
eight detector planes placed symmetrically around the Detector Solenoid axis, plus 



eight more planes projecting radially outward from each vertex of the octagon. 
Each plane consists of resistive straws approximately 300 cm long. A hit position 

in the radial and azimuthal direction is determined by the straw position and the 
drift time on the anode wire. The hit position in the axial direction is determined 
by the centroid of the imaged charge from the anode wire, as collected on cathode 
pads. 

The transverse tracker would have several advantages in comparison with a lon- 
gitudinal one since: 

• Normal non resistive straws are used without a pad system. 

• Gas manifolds, straw end-caps and chamber supports are all outside the conver- 
sion electron trajectories. 

• Shorter straws (0.7 - 1.3 m) are technically easier to build and are more robust 
against instability than longitudinal straws of 3 m. 

• Complications from the small tilt of each plane of the longitudinal tracker with 
respect to the magnetic field are avoided. 

• Transverse geometry provides a simple signature of an event since charged parti- 
cles cross a single straw only once. 

• There is a significant simplification dealing only with single chamber hits points. 

• The average number of spatially separated hits is a few times greater in com- 
parison with the number of spatially separated clusters in the longitudinal tracker. 
The hits are distributed uniformly in lobes along the transverse tracker. 

• The transverse tracker presents less material in the case of 15 yum straw wall 
thickness. The effective total thickness is 29*(15+2)*3.14 = 1550 /xm versus longi- 
tudinal tracker 8*3.5*25*3.14 = 2200 //m. 

• The pattern recognition can be performed with good precision even without drift 
time and amplitude information. 

• The Cu layer covering straw tubes suppresses significantly gas diffusion through 
a straw wall. 

• Capacitive crosstalk between channels is small. 

• Low energy Compton electron backgrounds have distinguished signature. 



Ill BACKGROUND SIMULATION 

Our analysis is based on a full GEANT simulation taking into account an indi- 
vidual straw structure. Multiple scattering and energy loss are taken into account. 
Isobutane {C4H10) gas is assumed to fill the tubes. 

For the subsequent analysis we take only events satisfying the following criteria: 
number of hits in the tracker is greater than 15; MC simulated energy release in the 
calorimeter is greater than 80 MeV; pitch angle is greater than 45 That leaves 
about 35% of the original events. 

The primary sources of charged particles in the tracker detector during the de- 
tection time are protons, neutrons and photons from muon capture by '^^Al nuclei 
and electrons from muon decay in orbit. The average and peak tracker rates from 
different backgrounds are presented in Table 1. The main source of background 
to muon conversion is muon decay in orbit (DIO). 

To study the tracker performance in the presence of the background the number 
of muon conversion events (10^) with initial momentum of conversion electrons 

A 



TABLE 1. The Average and Peak (in parenthesis) Tracker Rates 



Processes 


proton 


neutron 


7 


DIO 

< 55 MeV 


DIO 

> 55MeV 


Particles / process 
Particles / nsec 


0.1 

120(170) X 0.1 


1.2 

120(170) X 1.2 


1.8 

120(170) X 1.8 


0.9945 
80 (115) 


5.5-10-'' 
80 (115) 


Prob. Particle to Hit Tracker 


1.08-10-^ 


0.92- 10-^ 


0.9 • 10-^ 


4.0- IQ-" 


1.4- 10-^ 


Events in Tracker/30ns 


3.9 (5.5) 


4.0 (5.6) 


5.9 (8.2) 


1.0 (1.4) 


0.2 (0.3) 


Straws on in Tracker/event 


14.0 


5.2 


31.0 


15.7 


4.3 


Straw Rate (kHz) 


140 (200) 


53 (75) 


470 (650) 


40 (56) 


2.2 (3.3) 



produced in the target (105 MeV/c) was simulated and saved to a data file. The 
expected number of DIO events during the experiment time (10^ sec) in this region 
is 5.2 ■ 10^ events above 100 MeV/c. 

To study the main background the number of simulated DIO events in this energy 
region was chosen to be ten times greater. Five different random backgrounds 
(protons, neutrons, photons and DIO (see Table 1) were generated and saved to 
data files. 

The numbers in the table are the result of a GEANT simulations at the antici- 
pated beam intensity, with no artificial increase as a safety factor. No suppression 
was added for rejection of heavy ionization particles. 

Background was taken from the data files and added to conversion events and 
DIO events above 100 MeV/c. For example, we simulate the number of events N for 
the proton background according to a Poisson distribution with average 3.9 taken 
from background Table 1. Then we randomly pick-up N accidental proton events 
from the corresponding data file and add these events to the muon conversion or 
DIO event. In the same way the above procedure is repeated for all background 
types. 

In Figure 2 distributions in the number of real and background tracker hits are 
shown. The average number of real hits is 29. The real hits distribution starts 
from 15 hits, corresponding to the cut-off in the minimum number of selected hits. 
The maximum number of real tracker hits can reach 60. The average number of 
background hits is 300 which corresponds to an average straw rate about 800 kHz. 
The distribution in the number of background hits is a broad one and the maximum 
number of hits is about 900. 
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FIGURE 2. Distributions in the number of real and background tracker hits. 



IV PATTERN RECOGNITION 



There are twelve views in the tracker, each separated from the previous by 30 
degrees. In a typical event, a two-dimensional projection of the helical trajectory, 
a sine curve, is observed in three or four of the views. There are approximately 10 
hits in each view, and 29 hits in the event. The hits are grouped in lobes (see Figure 
3) with a typical gap between lobes about 60 cm (in these gaps an electron travels 
in vacuum). Therefore in the sensitive area of the tracker one gets track segments 
but not a complete track and the problem of combining information from different 
lobes arises. Since in a given view at least 4 hits are required to reconstruct four 
parameters describing a two-dimensional projection of helix often it is not possible 
to get helix parameters for a track segment in a single lobe and than to apply track 
element merging strategy. Quite opposite one has to use hits from different lobes 
for pattern recognition and then for momentum reconstruction. We will show that 
an appropriate mathematical approach allows to resolve the problem.. 

The sinusoidal projection of helix is described by four parameters: x'q , z'q, 
Rt, 

x[=x', + RtcosC-^) (1) 

where importantly, Rl and Rt-, related to the longitudinal and transversal mo- 
menta, and z coordinate are common to all projections. Coordinate x' is defined 
in a system for each given view. 

A simple and powerful reconstruction strategy starts from a single chamber view. 
To determine the parameters in this view a minimum of four hits are required. The 
total momentum and momentum pitch angle are determined. 




FIGURE 3. Lobes in the tracker. 



To find the five parameters describing the full 3-dimensional helix, one additional 
hit outside the plane is required: 



x[ = x'o + Rt cos(^^); y[ = y', + RT M^^) (2) 



The mathematical procedure to find the parameters is the following one: 



• four hits from a single view are taken to give the system of equations for hehx 
parameters in the coordinate system related with the view. 

• the system of four equations is reduced to a unique equation for Rl. 



Ax24'S'43'S'4i5'3i + Axi45'435'425'23 + Ax345'4i5'425'i2 = (3) 

where Axij — x'^ — Xj , Sij = sm{{zi — Zj)/2RL). 

• the equation is solved for numerically by the Newton-Raphson method. 

• the parameters are expressed in terms oi Rl. 

/ / / p V _ Aa:i45'43 sin((z4 + Z3)/2Rl) - Aa:345'4i sin((2:4 + zi)/2Rl) 
anl^o/ L) - ^^^^^^^ ^^^^^^^ ^ z^)/2Rl) - AX34-S41 cos((^4 + z^)/2Rl) ^ ' 

R, 



Cl - C2 
/ _ — X1C2 

Xq — 

Cl - C2 

where c = cos((2;i — Zq)/Rl). 

• a fifth hit from a different view is included to obtain the remaining parameters 
of helix. 



Xo — — 7 r- < xL sin a" — Ix'l — Rt cos( — a.")] sinfo;') > (5) 

sm(Q;" — a') I Rl /j v /j 



2 ( Zq ^ 

^0 = ■ / „ ;t \ x'q cos a" - [x^ - Rt cos( ^ ° - a")] cos(a') \ . (6) 

sm(Q;" — a') I Rl J 

where a' {a") is the angle of rotation from the global system to a given chamber 
system, 

zo^ Zq - a'RL- (7) 

It is important to note that the consideration of several lobes (see Figure 3) 
allows calculation of the longitudinal radius Rl with high precision AP/P Ri 10"^ 
due to the large distance between the lobes in comparison with their size. 

A two stage procedure was developed to provide the pattern recognition in the 
transverse tracker: 

• pattern recognition without drift time. At this step only information on 
centers of straw hits is used. The parameters of the reconstructed average helix are 
considered as a starting point to find an approximate helix which fits the straw hit 
centers. 

• pattern recognition with drift time. A deterministic annealing filter DAF 
[8] is applied to make a final background suppression and provide a starting point 
for track momentum reconstruction by the Kalman filter technique [9], [10]. 
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Pattern recognition without drift time information 



In the first level of the analysis procedure (the straw hit center approximation) 
only centers of straw hits were used as chamber hit coordinates x', z. It was assumed 
that there is only one useful muon conversion or DIO track. 

It is important to note that the straw hits allow us to reconstruct the helix 
parameters with high precision even without drift time information. Indeed, the 
uncertainty in the helix radius ARt is approximately a ~ D/y/T2 ^ 0.14 cm for 
the diameter of a straw tube D =5 mm. The average helix radius for muon con- 
version events is Rt ~ 25 cm . We can therefore expect the momentum resolution 
AP/P ~ 0.5 — 1% for a straw diameter of 5 mm. It is worth noting that using 
straw tube center positions without drift time information we can reconstruct the 
total momentum with the standard deviation a = 0.45 MeV/c for a straw diameter 
of 5 mm. 
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FIGURE 4. Momentum distribution of electrons entering the tracker. 



The initial momentum distribution of electrons entering the tracker and satisfying 
the criteria mentioned above is shown in Figure 4. This distribution is the a narrow 
one and can be fitted by a Gaussian with a standard deviation a = 0.21 MeV/c and 
average momentum Pj„ = 104.4 MeV/c in the range 104.2 MeV/c - 105 MeV/c. 
The width of the initial distribution of electrons is small in comparison with the 
total width of the momentum distribution reconstructed in the straw hit center 
approximation. 

In Figure 5 an initial map of real and background tracker hits in all tracker's 
views is shown as an example for a first simulated event. We will call it the sample 
event in the following. The number of real tracker hits is 29 and the number of 
background hits is 260 in this event. The real hits are depicted as cyan spots and 
the background ones are shown as red spots. Note that in the figure all 12 views 
were joined in one for all 18 tracker modules. 

One can see from the initial distribution of tracker hits that some of background 
hits are concentrated along straight lines ( Figure 5) providing a typical signature 
of Compton electrons from photon background. 
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FIGURE 5. Plot of real + background tracker hits before the selection procedure. 

The procedure of the pattern recognition and fitting without drift time in the 
presence of background can be described by the following steps: 

1) Rejection of hits produced by Compton electrons. 

• A grouping containing 6 modules is taken for a given view. If in this grouping 
there are A^^i > 2 hits with the same tube number (same straw (x') coordinate) the 
corresponding hit numbers are recorded. If, in a second grouping displaced by one 
module from the first, there are N2 > 2 hits with the same tube numbers as at the 
initial position of the grouping then hits from the first and the second grouping 
are marked as Compton electron hits. The procedure is repeated in forward and 
backward directions. At the end of this step all marked hits are temporary removed 
from the following analysis. This selection typically reduces twofold the number of 
background hits. Also about 5% of the real hits are rejected by the procedure but 
later they can be restored. 

Figure 6 demonstrates the remaining tracker hits after application of the rejec- 
tion procedure to the sample event. The suppression of background hits is evident. 
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FIGURE 6. Plot of real + background tracker hits after the first selection procedure. 



2) Selection based on calorimeter and tracker hits. 



To reconstruct the helix parameters four hits are needed. To reduce the num- 
ber of four hit combinations the calorimeter hit is used in each combination. The 

calorimeter coordinate resolution is a = 1.5 cm which is ten times greater than 
tracker hit resolution in the straw hit center approximation. In spite of the poor 
coordinate calorimeter resolution it is possible to get a significant background sup- 
pression. 

This procedure is described by the following steps: 

• For a given view all possible three hit combinations are chosen and four hits 
are formed by adding a calorimeter hit to the three hit combinations. The mathe- 
matical approach described above is applied to four hits. The algorithm allows us 
to calculate R , cos 9 , z'q , x'^m the chamber system of coordinates. Only the com- 
binations that survive a cut-off in the particle momentum pmin = 94MeK/ c < p < 
Pmax = ll4:MeV/c and pitch angle 6'inin — 41.4° < 9 < 6'inax = 66.4° are retained 
for subsequent analysis. 

• A fifth hit from a different view than the four hit combination is added. This 
allows us to calculate all 5 helix parameters in the global system of coordinates. A 
cut-off in xq and yo is applied at this step to five hits. The coordinates Xo,yo are 
required to be in an acceptable range —40 < Xo,yo < 40 cm, where the range is 
determined by the GEANT simulation 

• We look for hits correlated in 3 dimensions to choose good five hit combina- 
tions. For each selected five hit combination using the found helix parameters we 

reconstruct the helix, calculate all crossings of the tracker chambers for the recon- 
structed trajectory and define the total number of crossings N. Then we calculate 
how many times M at least one tracker hit matches the crossing in a road ± 3 cm. 
In an ideal case M should be equal to N. 

• The probability is evaluated to get M hits in the road. The probability that M 
of N crossings are in the road is estimated by a trial function 

Prob = — e*^(i - e)^-^ 

M\{N-My. ^ ' 

where e is the probability of that there is a hit in the road within ±3 cm. 

If the total probability Prob is greater than the threshold probability then the 
given helix is considered as a good candidate. The threshold probability and e 
were found empirically to be 0.001 and 0.95, respectively. 

As an output of the previous steps we get a collection of valid five hit combina- 
tions (helices). An individual hit is kept if it is in any valid hit combination, pro- 
viding a list of good tracker hits. It is worth emphasizing that, due to the strong 
spatial correlations between the helix hits in comparison with the un-correlated 
background, the number of false hits is reduced drastically on applying the road 
requirement. On average the number of background hits is reduced due to this step 
by a factor 20. 

In Figure 7 the surviving hits are shown. The significant reduction of background 
hits can be seen by comparing Figure 6 and Figure 7. 

Note that at this step no real hits are removed for the sample event. 

3) Selection based on tracker hits. 

1 n 
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FIGURE 7. Plot of real + background tracker hits after the second selection procedure. 



In this selection procedure the tracker hits found in step 2 are used. This proce- 
dure repeats the previous one with the following changes: 

• Hits are chosen only from the tracker hits which allows to improve the perfor- 
mance of the selection procedure since tracker hit coordinates are defined signifi- 
cantly better than calorimeter hit coordinates. 

A plot of the number of four hit and five hit combinations in the presence of the 
background (in average about 120 and 630 respectively) is shown in Figure 8 (a) 
and (b). 
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FIGURE 8. Number of four hit and five hit combinations. 



As seen in Figure 8 these distributions are broad ones reaching a few thousand 
combinations. 

• For each five hit combination selected above we evaluate the position of a hit 
in the calorimeter Xevai , T'eval on the basis of the defined helix parameters. The 
hit combination is accepted only if the evaluated hit matches the calorimeter hit 
within a road : \xcaio — Xevai\ < Tcm . 

• Parameters at this step are more restrictive: the road width is taken to be 0.75 
cm and the probability to be within the road e = 0.97. 

Figure 9 shows tracker hits remaining after this step for the sample event. The 
number of background hits is reduced typically by a factor 4-5 due to this step. 

4) Selection including restored tracker hits. 

1 1 
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FIGURE 9. Plot of real + background tracker hits after the third selection procedure. 



• Some real tracker hits are lost at previous steps of the procedure especially 
at Compton electron hits rejection. To restore the lost hits an average helix for a 
given event is reconstructed and tracker hits that match the average helix in a road 
±1 straw are added to the list of valid tracker hits obtained above. 

• Step 3 of the procedure is repeated for the extended list of tracker hits. On 
average the list of valid hits is extended by 1 real and 2 background hits. 

• The union of all hits in valid five hit combinations in the extended list is used 
to provide an input for an average helix. A global fit is applied to reconstruct 
the helix parameters on the basis of the list of the selected tracker hits and the 
parameters of the average helix are considered as a starting point for the fit. 

In this section the pattern recognition procedure without drift time based on 
the straw hit center approximation was developed. In this approximation no am- 
biguity due to mirror hits arises and the total momentum is reconstructed with 
the standard deviation a = 0.45 MeV/c. The selection procedure based on 3D 
space correlations between real tracker hits inside the road significantly reduces 
background. The overall background rejection factor is about 130 for the pattern 
recognition procedure without drift time. 



Pattern recognition with drift time 

The second stage in the pattern recognition procedure uses the fitted helix and 
drift time for hits selected at the first stage. The reconstruction of helix parameters 
can be improved by taking into account that in addition to a straw coordinate for 
each hit a chamber gives the radius r calculated from the measured drift time 
fneas ^ The crrors (a) in radius measurements were taken to be 0.2 mm. This 
radius r carries an ambiguity as to whether the track passed left or right of the 
wire. The search of two possible up and down hit positions lying on the circle of 
radius r is based on the fitted helix obtained previously in the straw hit center 
approximation. We can call these up and down points as true and mirror ones. Up 
and down points are extracted (see Figure 10) from the intersections of a normal 
to the helix through the straw center and the circle of the drift radius r. In this 
case coordinates of up and down points are given by 

1 o 



X, = XcT rsinP; z = Zc±rcosP; tan(3 = —{Rt I RL)sin{{zc — z'q) / Rl) (8) 
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FIGURE 10. Reconstruction of up and down points by the helix obtained in the straw hit 
center approximation. 

The uncertainties As (see Figure 11) in the determination of up and down point 
positions are small and can be evaluated in the following way. The direction of 
2D helix in the chamber coordinate system is given by tanj3 (see Figure 11 and 
Eq.(8)). 
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FIGURE 11. Uncertainties in up and down point positions. 

Since the uncertainty in Rt is the dominant one the uncertainty in the direction 
of the helix is of the order A/5//? ~ ARt/Rt ~ AR/R ~ AP/P where AP/P < 
0.01, /3~ 1. 

The uncertainty in the determination of up or down point position 
As ~ A/? r < j3{AP/ P)rtube ~ 1 x 10~^ x 0.25 cm ~ 25/im which is much less than 
a measurement precision 200 /xm. 

To reject background hits remaining after pre-selection stage and to resolve the 
up - down ambiguity the deterministic annealing filter (DAF) and the Kalman filter 
(KF) are applied. 

In principle at this stage the Kalman filter (KF) approach [10] could be applied. 
However an application of the Kalman filter (see Appendix A) requires that the 
problem of assignment of hits to a track has been entirely resolved by the preceding 
selection procedure. If this is not the case the filter has to run on every possible 
assignment choosing the best one according to the chi-square criterion. For the 
number of tracker hits greater than 15 this combinatorial search is computation- 
ally expensive and practically unfeasible. Therefore as the last step of pattern 
recognition we will use the deterministic annealing filter (DAF) [8]. 

1 Q 



DAF is a Kalman filter with re- weighted observations (see Appendix B) . For the 
DAF procedure we introduce artificial layers placed at the chamber straw centers 

in order to have competing true and mirror points in one layer. To overcome the 
problem of insufficient information in the initial phase of the filter, an iterative 
procedure is applied. After a first pass of filter and after smoothing, the track 
position can be predicted in every layer of the tracker. Based on these predictions, 
the assignment probabilities for all competing hits can be calculated in every layer. 
If the probability falls below a certain threshold, the hit is excluded from the 
following consideration. The assignment probabilities of the remaining hits are 
normalized to one and used as the weights in the next iterations of the filter. 
In our case the operation of DAF is described by the following steps: 

• every true and mirror point is projected on a layer corresponding to the center 
of straw in the direction defined by the fitted helix for a given event . 

• initial probabilities of competing points in the layer are assumed to be equal. 

• annealing schedule is chosen according to the following formula 

Vn = V {j^ + 1) ior a, variance of observations, where V = cr^ and a — 200 fim. The 
annealing factor f is chosen to be either 1.4 or 2. 

• standard Kalman filter runs on all layers taking observations as weighted mean 
according to assignment probabilities. 

• the filter runs in the opposite direction, using the same weighted mean as the 

forward filter. By taking a weighted mean of the predictions of both filters at every 
layer, a smoothed state vector and its covariance matrix are obtained. 

• based on these predictions and the covariance matrix, the assignment proba- 
bilities of the hits are calculated. If combined hit probability for true and mirror 
point falls below a certain threshold (10^^), the hit is rejected. The assignment 
probabilities of the remaining points in the layer are normalized to one and used 
as the weights in the next iterations of the filter. 

• iterations in n stop if one of the following conditions is satisfied: 

1) for the reconstructed track > xLax where xLax = 1000 or < 0. 

2) for the reconstructed initial momentum Pj„ < 94MeV"/c or P^n > llAMeV/c 

3) variation of is small in comparison with the previous iteration \Xn+i~Xn\/Xn < 
0.01. 

• when iterations stop (on average after 7 iterations) The DAF procedure is 
repeated for the different annealing factor f. From results corresponding to two 
annealing factors we choose that one corresponding to the minimum x^- 

Figure 12 shows the distribution of tracker hits for the sample event (recall 
that initially we had 260 background hits for this event). There are no surviving 
background hits and missing real hits in this case. Reconstructed lobes are clearly 
seen in the Figure. We conclude that DAF is effective in rejection of background 
hits remaining after pre-selection procedure and also it provides a good starting 
point for the track reconstruction. 

1 A 



2D INITIAL REAL S BACKGROUND TRACKER HITS 



900 950 1000 1050 1100 1150 



FIGURE 12. Plot of real + background tracker hits after the DAF selection procedure. 



Figure 13(a) represents the number of real hits lost by the pattern recognition 
procedure. Some of the real tracker hits (0.8 hits ~ 2.7%) are lost due to the 
selection procedure. 

Figure 13(b) represents the number of background hits remaining after the 
selection procedure. As one can see from Figure 13(b) the number of background 
hits remaining is 0.38 hits in comparison with the primary 300 hits. So a total 
background suppression factor is 300/0.38 ~ 800. 
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FIGURE 13. Distributions in the number of missing real hits and remaining background hits. 



V RECONSTRUCTION BASED ON KALMAN FILTER 

Background studies presented in this subsection are based on the application of 
the Kalman filter to hits selected by previous pattern recognition procedure. 

At this stage the reconstruction is based on the hits selected by the pattern 
recognition procedure described above. 

In principle we could use the results obtained by DAF for the track reconstruction 
but our analysis shows that an application of a combinatorial Kalman filter to hits 
selected by DAF provides better precision. 

In the last decade the Kalman filter (KF) approach [10] has been extensively 
exploited for track fitting in high energy physics. This approach possesses the fol- 
lowing features for effective track fitting: 

• multiple scattering and energy losses are included in a natural way; 
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• a 3D trajectory is restored that approximates closely the real one; 

• complex tracker geometries are handled in a simple way; 

• NxN matrix inversion, where N is the total number of measurements, is avoided; 

• control for error propagation is provided; 

• trajectory is reconstructed progressively from one measurement to the next, im- 
proving the precision with each step; 

• initial and final momenta of a particle crossing the tracker are reconstructed. 

The KF is very useful because it simultaneously finds and fits the track; it is 
much more economical than the conventional least-squares global fit. The KF 
is a "progressive" step by step method whose predictions are rather poor at the 
beginning of the track at the first stage of filtering. Since a state vector = 
{x,y,tx,ty,l/pL) at point k (see definitions in Appendix A) has five parameters 
we need approximately ~ 6-7 straw hits to get good KF prediction precision. The 
prediction step, in which an estimate is made for the next measurement from the 
current knowledge of the state vector, is very useful to discard noise signal and hits 
from other tracks. Assuming the validity of the helix track model for each step, 
the KF propagates the track in 3D space, from one 2D surface to the next. 

Below we will use the standard notations: 

x}^_,_-^ is a prediction, i.e. the estimation of the "future" state vector at position 
"k+1" using all the "past" measurements up to and including "k". 

x}^ is a filtered state vector, i.e. the estimation of the state vector at position "k" 
based upon all '"past" and "present" measurements up to and including "k" . 

The same notations will be held for the covariance matrix C, noise matrix Q and 
so on. 

The Kalman filter algorithm can be divided into three major steps. 
INITIALIZATION 

The KF can start from arbitrary parameters and an infinite covariance matrix 
but for track finding applications it is significantly better to fix somehow the initial 

state if possible. 

The initialization of forward and backward KF algorithms is quite simple. It 
starts from the artificial point and initial parameters calculated from the previ- 
ous stage of reconstruction procedure. The initial covariance matrix is empirically 
found to be diagonal with matrix elements being much greater than the correspond- 
ing uncertainties: = (0.3,0.3,0.03,0.03,0.000003). 

FILTERING 

Once the KF is initialized it makes standard consequent steps. The current state 
at the k-th step is defined by the state vector Xk, the state covariance matrix 
and the current straw hit. To take the (k+l)-th step it is necessary to: 

• update x{^, to take into account the ionization losses 

• define the next hit object (up or down point) 

• propagate the parameters and the covariance matrix to: x}^ C}^ — > 

• update Ck_,_i to take into account multiple scattering C^_^_i C^+i + Qk 

• calculate the Kalman matrix Kk+i 

• update the covariance matrix C^+i 

• calculate residuals r^'^\ and their covariance matrices R^^i 




• store all information defining the new state 



SMOOTHING 

In the standard Kalman filter algorithm the smoothing is a well defined proce- 
dure. Smoothing allows one to obtain the best estimate of the track parameters at 
any trajectory point using all hit information accumulated during the KF propa- 
gation. 

Figure 14 (a) displays the results of the KF forward filtering for the total momen- 
tum Ptot reconstruction at each tracker hit position. Figure 14 (b) displays how 
the KF smoother, based on all hit information accumulated during the KF filtering, 
improves the total momentum Ptot reconstruction at each tracker hit position. 
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FIGURE 14. The total momentum reconstructed by the forward Kalman filtering a) and by 
the smoothing b) for each hit of a selected event. The y-axis is Ptot, the x-axis is the z-ordered 
hit number. 

Figure 15 (a) displays the results of the KF backward filtering for the total 
momentum Ptot reconstruction at each tracker hit position. Figure 15 (b) displays 
how the KF smoother, based on all hit information accumulated during the KF 
filtering, improves the total momentum Ptot reconstruction at each tracker hit 
position. 
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FIGURE 15. The total momentum reconstructed by the backward Kalman filtering a) and by 
the smoothing b) for each hit of a selected event. The y-axis is Ptot, the x-axis is the z-ordered 
hit number. 

The Kalman filter approach is effective in the resolving of up - down ambiguity. 
By applying the Kalman filter at each step up and down points are considered and 
the point providing the best for the trajectory is selected as the true point. This 
procedure is approximately linear in the number of tracker hits in comparison with 
the a combinatorial search, which is not feasible. 
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At this reconstruction stage due to left-right ambiguity for the straw drift cham- 
ber we have a set of true and mirror points for each straw hit. As discussed above 

the reconstruction procedure for each tracker hit defines true and mirror points 
and one of them is close to the real point with high precision ~ 25 fim. So we 
can formulate our goal as to find a true point combination for N straw hits in the 
presence of N mirror points. 

The procedure to find the best approximation to the true point combination for 
N straw hits in the presence of N mirror points for one track from a muon-electron 
conversion event is based on the following steps: 

• chose the first eight straw hits and built 2^ (256) possible hit combinations 

corresponding to up and down points. The KF forward and backward procedures 
described above are applied to these combinations. Only those combinations which 
satisfy a rather loose cut, < 30, are retained (typically about 10 out of the 
initial 256 combinations); 

• make a loop for all retained combinations with fixed up and down points for 
the first eight straw hits. For the 9th and higher straw hit, the up and down point 
choices are take into account and the point with minimal incremental for this 
point is selected for the further KF propagation step. If both incremental satisfy 
the cut < 10 the second point is stored in the stack to make an iterative loop. 
A single combination of all possible steps defines a candidate track. At this stage 
on average 45 candidate tracks are stored in the stack per event; 

• make a loop for all combinations from the stack. For each new hit added to the 
hits restored from the stack again the two up and down point choices are taken into 
account and the point with the minimal incremental foi' this point is selected 
for the further KF propagation step; 

• select up and down track combinations with the minimal for the track; 

• select a track satisfying the cut < 70; 

• select a track with the difference between the forward (Pin_f) and backward 
(Pin_b) reconstructed input momentum satisfying the cut |Pm_/ — PinJb\ < 0.7 
MeV/c. 

The Kalman filter reconstructs a trajectory of a particle in three dimensions. The 

trajectory is bent each time it crosses a tracker plane due to multiple scattering. 
Therefore, the reconstructed track is a set of helices that intersect at the planes. 
This is the track followed by the particle. 

Figure 16 displays a 2D projection of 3D trajectory reconstructed by the Kalman 
filter for the sample event. As above in this figure all 12 views were joined in one 
for all 18 tracker modules. The 2D trajectory is shown only for sensitive area of the 
tracker and for each view the trajectory is in a different color. For the sample event 
real hits are in four tracker's views and the reconstructed lobes for these views are 
clearly seen in the figure. 

Due to a scale in this figure the 2D trajectory looks as an ideal sine curve and 
tracker hits look like spots of different size. In order to see a detailed behavior 
of the trajectory and hit positions a dynamical zoom is applied to a rectangular 
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region indicated in Figure 16 in the x-range 49-50 cm and z-range 1029-1029.5 cm. 

Figure 17 demonstrates the magnified region of the tracker. Blue fine in the 
figure represents the reconstructed trajectory. Two circles represent two hits in 
chamber. These circles look like ellipses due to different axis scales. 
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FIGURE 16. Plot of real + background tracker hits after the DAF selection procedure. 
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FIGURE 17. Zoom enlargement of the region in Figure 16. 



Radii of the circles are directly proportional to the drift times. For illustrative 
purposes two nearest to the straw center points corresponding to the real trajectory 
obtained in Monte Carlo simulation are shown in the Figure in the form of diamonds 
(note that in the pattern recognition procedure only radii were used but not these 
points). The measurement uncertainty was taken into account assuming that circle 
radii are distributed normally about the simulated radii with a = 200yum. For this 
reason the position of one of the nearest points is not on a circle in Figure 17. 

In Figure 17 due to the corresponding scale the trajectory looks like a straight 
line. It is tangent to one of two circles obtained on the basis of drift time. In the 
region under consideration the deviation of the trajectory from the nearest point 
is less than 0.2 mm. The change in the direction of the trajectory due to multiple 
scattering can not be seen in the Figure because of the smallness of the average 
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angle of the scattering. 

Figure 18 shows transverse xy-projection of the trajectory for the sample event. 
In this projection the trajectory looks approximately as a circle. 
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FIGURE 18. The transverse projection of the 3D trajectory reconstructed for the sample event. 

However if a specific region of the tracker is magnified by the dynamical zoom 
one can see in Figure 19 that the shape of the circle is distorted due to multiple 
scattering and energy loss. More than two turns of trajectory are clearly seen in 
the Figure. 
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FIGURE 19. Zoom enlargement of the region in Figure 18. 

Note that by using the KF filter a momentum of a particle can be reconstructed 
at any point of the tracker. For our purposes the most important is the momentum 
of a particle entering the tracker, which in the following we will call the input 
momentum. 

Figure 20 demonstrate 3D trajectory reconstructed for the sample event. The 
trajectory looks in this scale as a helix, but we remind that it consists of many 
helix parts. Also in the Figure tracker's hits generated by Monte Carlo simulation 
program are shown. 

The distribution in the difference between the initial momentum (Pin_f) recon- 
structed by the Kalman filter and the generated initial momentum (Pin) is shown 
in Figure 21 in linear (a) and logarithmic (b) scale. According to this distribution 

on 



3D HELIX RECONSTRUCTED 



BY KF & 3D TRACKER HITS 




1200 



800 



FIGURE 20. 3D trajectory reconstructed for the sample event. 

the intrinsic tracker resolution is a = 0.12 MeV/c if one fits the distribution by a 
Gaussian in the range -0.3 - 0.7 MeV/c. 
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FIGURE 21. Distribution in the difference between the input reconstructed momentum based 
on the Kalman filter and the simulated input momentum with background. 

Figure 22 shows a distribution in the input momentum (Pin_f) reconstructed 
by the Kalman filter in linear (a) and logarithmic (b) scale with background. This 
distribution is characterized by the standard deviation a = 0.25 MeV/c of the 
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FIGURE 22. Distribution in the input momentum reconstructed by the Kalman filter with 
background. 

reconstructed input momentum for a Gaussian fit in tlie range 104-106 MeV/c. 

Note that the trajectory reconstructed by the Kalman filter consists of many 
helix parts. The reconstructed input momentum resolution by the Kalman filter is 
a = 0.25 MeV/c. This resolution significantly better than the resolution a = 0.35 
MeV/c obtained by a single helix fit. 

The overall reconstruction acceptance is 22.1 % for muon conversion events with 
the momentum above a threshold momentum of 103.6 MeV/c. 

Comparing these results with the results of the reconstruction without back- 
ground we get the difference in tracker resolution 1.5 % and the difference 2.7 
% in overall acceptance (see Appendix C). Therefore the tracker resolution and 
overall acceptance are not affected significantly at the considered background level 
(explained in Section 3). 

A summary of the critical selection criteria used in the momentum reconstruc- 
tion is shown in Table 2. The efficiencies are for the selection criteria applied in 
consecutive order. 

An overall acceptance for muon conversion events with momentum above thresh- 
old momentum Pth is 22.1% . We define a threshold momentum, above which events 
are considered as the useful ones by Pth = Pmax — where Pmax = 104.3 MeV/c 
the most probable reconstructed momentum. If A = 0.7 MeV/c is chosen then Pfh 
= 103.6 MeV/c. 



TABLE 2. A summary of the reconstruction se- 
lection criteria 



Selection criterion 


Efficiency 


Calorimeter energy above 80 MeV 


0.53 


Required pitch angle at the tracker 


0.86 


At least 15 hits in the tracker 


0.87 


Position match in the calorimeter 


0.95 


Requirements on fit quality 


0.79 


Detected energy above 103.6 MeV 


0.74 


Overall acceptance 


0.22 
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The main factors entering into the experimental sensitivity are the running time, 
the proton intensity, the probabihty per proton that a is produced, transported 
and stopped in the stopping target, the fraction of stopped muons that are cap- 
tured (as opposed to decay), the trigger efficiency and the tracker reconstruction 
acceptance. We do not include in this table loss of events due to accidental cosmic 
ray vetoes, dead-time losses and losses due to straw chamber inefficiencies, all of 
which are expected to be small. According to our analysis taking into account the 
straw efficiency 97% the overall acceptance is reduced from 22.1% to 21.7% . 

Table 3 shows expected MECO sensitivity for a one year (10^ s) run. 



TABLE 3. A summary of the expected MECO sensitivity. 



Running time (s) 
Proton flux (Hz) 


10^ 
4 • 10" 


Probability of /i/p transported and stopped in target 

capture probability 
Fraction of /i which are captured in time window 
Trigger efficiency and the selection criteria 


0.0025 
0.6 
0.49 
0.22 


Detected events for R^e ~ 10^^'' 


6.5 



Muon DIO events are the most important background for the experiment. The 
main background from muon DIO events in the presence of background tracker 
hits was simulated and reconstructed. Based on the simulated DIO events in the 
momentum range above 100 MeV/c the track pattern recognition and momentum 
reconstruction were performed in the presence of background tracker hits (protons, 
neutrons, photons, DIO) by applying the selection criteria discussed above. 
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FIGURE 23. Distributions of DIO reconstructed momentimr in the presence of the background 
above 100 MeV/c . 
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The number of primary DIO events simulated in the momentum range above 100 
MeV/c was 10 times more than expected and 3 background events were found (see 
Figure 23). So the background is expected to be 0.3 events. 

It is important to note that input reconstructed momenta for these three back- 
ground events are very close to simulated input momenta. The input simulated 
momenta 103.48, 103.6 and 103.74 MeV/c for these events have to be compared 
with input reconstructed momenta 103.6, 103.76 and 103.8 MeV/c, respectively. 

At the present level of pattern recognition and momentum reconstruction studies 
we can expect the background from DIO events in the range above 100 MeV/c ~0.3 
events compared to 6.5 signal events for R^e = 10~^^. 

VI CONCLUSION 

A study of the impact of background on the performance of the transverse tracker 
proposed for the MECO experiment is presented. Background from capture pro- 
tons, neutrons and photons, and from muon decay in orbit was generated using 
GEANT3. The effective average straw tube rate from these sources was 800 kHz 
at the proposed muon beam intensity of 2 x 10^^/x~/sec. 

A pattern recognition procedure based on a Kalman filter technique was devel- 
oped to suppress background and assign hits to tracks. In the first stage of this 
procedure, straw hit center coordinates, without drift time information, were used 
to reduce the background by a factor ~ 130. In the second stage, the full drift 
time information and a deterministic annealing filter were used to obtain an ad- 
ditional six-fold suppression. The total suppression of 800 reduces the number of 
background hits on average from an initial value of 300 to approximately 0.4 per 
event. About 0.8 hits of the 29 real hits typically recorded, or 2.7 %, are lost in 
the process. 

It was found that in the presence of background the resolution of the tracker is 
(7 = 0.12 MeV/c and the overall setup acceptance for muon conversion events with 
momentum above the threshold momentum 103.6 MeV/c is about 22 %. At the 
considered background level the tracker resolution and the overall acceptance are 
not affected significantly by presence of the background: the tracker resolution is 
changed by 1.5 % and the overall acceptance by 2.7 %. 

Additional constraints on the background, not considered here, may be imposed 
if the drift time measurement is supplemented by a measurement of the pulse am- 
plitude at the anode wire. A straightforward, crude measurement of the amplitude 
is sufficient to reduce significantly the background from heavily ionizing particles, 
i.e., the capture protons that comprise 30 % of the background hits in the above 
study. We estimate too that a significant background suppression, of 10-20, would 
be achieved if resistive anode wires were used in place of conducting wires to ob- 
tain a measurement of the hit position along the wire. This would improve the 
resolution as well. 

At the present level of pattern recognition and momentum reconstruction studies 
the background from DIO events above 100 MeV/c is about 0.3 events. This is to 
be compared to 6.5 signal events for R^e = 10~^^. 

The study carried out shows that the developed procedures of pattern recogni- 
tion and momentum reconstruction in the case of the transverse tracker provide a 
required precision for lepton number violation search at a sensitivity level about 
10-1^. 
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APPENDIX A. KALMAN FILTER AND ITS 
APPLICATION TO TRACK FITTING 



The Kalman filter is an algorithm that processes measurements to deduce an 
optimum estimate of the past, present, or future state of a dynamic system by 
using a time sequence of measurements of the system behavior, plus a statistical 
model that characterizes the system and measurements errors, plus initial condition 
information. 

The Kalman filter addresses the general problem of trying to estimate at different 
points (1 < A; < n) the state Xk of a discrete process that is governed by the linear 
stochastic difference equation 

Xk = Fk_iXk_i + Wk-i (9) 

with a measurement mk that is 

mk = Hxk + £k- (10) 

The system equation (9) is not deterministic since the track experiences stochas- 
tic processes such as multiple scattering, bremsstrahlung, etc. These processes are 
taken into account by the process noise Wk . £k represents the measurement noise. 
Wk and Sk are assumed to be independent of each other with zero expectation 
values: 

-E{wk} = 0, cov{wk} = Qk , I < k <n, 
-^'{^k} — 0, covjek} — Vk , 1 < A; < n, 

where Qk and Vk are process noise and measurement noise covariances, respec- 
tively. 

Eq. (9) in the absence of the last term is the standard equation of motion with a 
propagator Fk-i (transport matrix). Note that at the moment Fk-i is assumed to 
be constant. 

Regarding a track in space as a dynamic system the filtering technique is applied 
to the track fitting. For example, in the case of a particle moving in magnetic field 
this can be done naturally by identifying the state vector Xk of the dynamic system 
with a vector Xk = (x, y, tan^x, tan6'y, 1/pl) of 5 parameters uniquely describing 
the track at each point of the trajectory. The F matrix propagates the state vector 
on one plane to the state vector on the next plane combining position information 
with directional information. The transport matrix implicitly contains information 
about a gap between planes. 

In general the set of parameters Xk is not measured directly; only a function of 
Xk , Hxk is observed. For example, in the case of the transverse tracker one does 
not measure Xk but xf — x cos a + y sin a in the chamber coordinate system which 
corresponds to 

H = (cos a, sin a, 0,0, 0) (11) 
There are three types of operations to be performed in the analysis of a track. 

• Prediction is the estimation of the "future" state vector at position "k" using 
all the "past" measurements up to and including "k-1". xj^"""^ is a prediction 
(a priori state estimation). 
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• Filtering is the estimation of the state vector at position "k" based upon all 
"past" and "present" measurements up to and including "k" . is a filtered 
state vector (a posteriori state estimation). 

• Smoothing is the estimation of the "past" state vector at position "k" based 
on all "n" measurements taken up to the present time, xj^ is a smoothed state 
vector. 

The first step to estimate Xk is the prediction (time update): 

xl:-^ = Fk-ix|::} (12) 

Cl:-^ = Fk_iCi::iF^_i + Qk_i (13) 

where Eq.(12) projects the state ahead and Eq.(13) projects the error covariance 
ahead. 

The filtered estimate (measurement update) x|^ is calculated as a weighted 



mean of the prediction and the observation: 

Kk = C|:-^H^[HkC|:-^H^ + Vk]-^ (14) 

Xk = + Kk[mk - Hkxi:-!] (15) 

C|: = [I-KkHk]Cl:-^ (16) 



Eq.(14) computes the Kalman gain matrix defining the correction to the predicted 
state due to the current observation. Eq.(15) updates the prediction with the 
measurement and Eq.(16) updates the error covariance. The error covariance may 
be also expressed in a computationally superior form 

Cl: = [I - KkHkjCl^-Ml - KkHk]"^ + KkVkK^. (17) 

The filtering is a recursive operation. The prediction step and the filtering step 
are repeated for the next plane proceeding progressively from plane "1" to plane 
"n" . The state vector at the last filtered point contains always the full information 
from all points. 

At each step one can calculate the filtered residuals r}^ , the covariance matrix of 
the filtered residuals Rk and the filtered x^'- 

rj^ = mk - Hkxj^ 
= Vk - HkCkHk 

xi - (R!:)-^ri: 

where xl - distributed with dim(mk) degrees of freedom. The total of 
the track is given by the sum of the xt contributions for each plane. 

The system of equations defining the Kalman filter represents an asymptotically 
stable system, and therefore, the estimate of the state vector xj^ becomes indepen- 
dent on the starting point Xq , Cq as k is increased. 
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When the last plane (nth) is taken into account the Kalman filter performs the 
final step which is a smoothing. The filter runs backward in time updating all 
filtered state vectors on the basis of information from all n planes. The equations 
describing the smoothing are given by 

Ak = CkFk (Ck+i)~^ 

— + Ak(Xk_|_i — 

= mk - Hkxl^ 

Rk = Vk - HkCkHk 

Until now it was assumed that the problem of estimation of a discrete-time pro- 
cess is described by a linear stochastic differential equation. However for example in 
the presence of a magnetic field the track propagator F is non-linear. Let's assume 
that the process of a particle propagation is governed by the non-linear stochastic 
differential equation 

Xk = f{xk-i) + Wk-i (18) 

with a measurement m in the form Eq.(lO). f is a non-linear function. The 

Kalman filter can be applied to this system by linearizing the system for example 
about the estimated trajectory. If deviations between the estimated trajectory and 
the actual trajectory remain sufficiently small the linear approximation is valid. 
The non-linear equation (18) can be written down in the linearized form as 

Xk = f (xj^li) + F • (xk-i - x^Zl) + Wk-i (19) 

where as before Xk,mk are the actual state and measurement vectors, x|^ is a 
filtered estimate of the state at step k. F is Jacobian matrix 

Fy = afi(xj^:i)/9xj (20) 

Therefore the complete set of extended Kalman filter equations is given by 
Eqs.(13)-(16),(19) by using F in the form (20). 

In order to apply the extended Kalman filter to a track fitting for a particle 
moving in uniform magnetic field (the magnetic field is in z direction) one has to 
choose the state vector parameters, define the initial state vector and calculate 
the transport matrix F, the projection matrix H, and the noise matrix Q. As 
it was mentioned above in this case the state vector can be chosen in the form 
Xk — {x, y, tx, ty, 1/pl) where x, y are the track coordinates in the tracker system, 
tx = Px/pl, ty = Py/PL define the track direction. The projection matrix H is given 
by Eq.(ll). Due to multiple scattering the absolute value of electron momentum 
remains unaffected, while the direction is changed. This defiection can be described 



using two orthogonal scattering angles, which are also orthogonal to the particle 
momentum [11]. In terms of these variables the noise matrix is given by 
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For the variance of the multiple scattering angle the well-known expression is 
used 



< >= {13.6MeV/pf[l + 0.038 ln(t/X^j)]t/XR 



(21) 



where Xji is a radiation length, t is a distance traveled by the particle inside a 
scatterer. Energy losses are taken into account by 



p/ = p- < dE/dx > t. 



(22) 



APPENDIX B. DETERMINISTIC ANNEALING FILTER 

Track reconstruction in modern high energy physics experiments faces a signif- 
icant amount of noise hits in a detector. The track fit thus is confronted with 
several competing hits in detector's layers. The Kalman filter (see Appendix A) 
is now widely used for the reconstruction of the track parameters in high energy 
physics. However the application of the Kalman filter requires that the problem of 
assignment of the detector hits to track candidates has been entirely resolved by the 
preceding selection procedure. If this is not the case, the filter has to run on every 
possible assignment to select the best one by chi-square criterion. Obviously this 
approach is computationally expensive and practically unfeasible for a considerable 
amount of noise hits. For this reason the Deterministic Annealing Filter (DAF) 
was developed [8]. In DAF there is an additional validation feature eliminating 
hits which are not compatible with the predicted track position. 

The deterministic annealing filter itself is a Kalman filter with re-weighted obser- 
vations. The propagation part of DAF is identical to the standard Kalman filter. 

The filtered estimate (measurement update) at layer k is calculated as a 
weighted mean of the prediction xj^"-*^ and the observations mj^., i = 1, 2, ...nt: 

xi: = xr ' + Kk E pLK - Hkxi:-^] (23) 

i=l 

where is the assignment probability of observation mj^. Kk is the Kalman 
gain matrix which is given by 

Kk = [[Crr' + PkH^Vi:^Hk]-^H^V^i (24) 

where pk is the sum over all weights Pk, Hk is the measurement matrix, Vk is 

the variance of the observations. 

The covariance matrix Ck of the updated estimate x{^ is written as 

00 



Cj:=[[Cn-' + PkH^V,^Hk]-^ 



(25) 



After completion of the forward filter a backward filter runs in opposite direction, 
using the same weights as the forward filter. By taking a weighted mean of the 
filtered states of both filters at every layer a prediction for the state vector x{J* 
along with its covariance matrix Cg* is obtained, using all hits except the ones 
at layer k. (The asterisk indicates that the information from layer k is not used 
in this prediction.) Initially all assignment probabilities for the hits in each layer 
are set to be equal but based on the estimated state vector xJJ* and its covariance 
matrix, the assignment probabilities of all competing hits are then recalculated in 
the following way: 

pL ~ ¥^(1111,; Hkxr , Vk + HkCrH^) (26) 

where (/^(x; ji, V) is a multivariate Gaussian probability density with mean vector 
H and covariance matrix V. 

If the probability falls below a certain threshold, the hit is considered as the false 
one and is excluded from the list of the hits assigned to the track. 

However at this step we cannot be sure in calculated probabilities especially in the 
initial phase due to insufficient information for the filter. This problem is overcome 
by adopting a simulated annealing iterative procedure. This is an additional feature 
of DAF. 

The simulated annealing optimization algorithm is based on an analogy between 
the behavior of a material heated past its melting point that is slowly cooled (an- 
nealed) to form a single crystal. If the cooling proceeds slowly enough, the crys- 
talline state reached at zero temperature will have all the atoms fixed in a perfect 
lattice structure, corresponding to the lowest possible energy of the system (global 
minimum) . 

In the same way in track fitting the simulated annealing allows to avoid a local 
minimum and find the global one corresponding to the minimum chi-square for the 
track. 

DAF annealing algorithm can be described in the following way. The annealing 
schedule is chosen for example in the form Vn = V(A//^ + l) where the annealing 
factor / > 1 and factor A» 1. This provides that the initial variance is well above 
the nominal value V of the observation error but the final one tends to V. After 
each iteration the assignment probabilities exceeding the threshold are normalized 
to 1 and used again as weights in the next iteration, and so on. The iterations 
generally are stopped if the relative change in chi-square is less than correspondent 
control parameter (typically of the order 0.01). 

Since we deal with the stochastic process the best result can be reached repeating 
the DAF procedure for a few different annealing factors f and then choosing the 
result corresponding to the minimum chi-square. 



APPENDIX C. TRACKER RESOLUTION 

This appendix demonstrates the results of apphcation of the pattern recognition 
and reconstruction procedure for the conversion events without background. 

The distribution in the difference between the input reconstructed momentum 
(Pin_f) based on the Kalman fiher and the simulated input momentum (Pin) is 
shown in Figure 24 in linear (a) and logarithmic (b) scale. According to this dis- 
tribution the intrinsic tracker resolution is a = 0.12 MeV if one fits the distribution 
by a Gaussian in the range -0.3 - 0.7 MeV. 




FIGURE 24. Distribution in the difference between the input reconstructed momentum based 
on the Kalman filter and the simulated input momentum without background. 
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FIGURE 25. Distribution in the input momentum reconstructed by the Kalman filter without 
background. 



Figure 25 shows a distribution in the input momentum (Pin_f) reconstructed 
by the Kalman filter in linear (a) and logarithmic (b) scale without background. 
This distribution is characterized by the standard deviation a = 0.25 MeV of the 



Q1 



reconstructed input momentum for a Gaussian fit in the range 104-106 MeV. The 
overall reconstruction acceptance is 22.7 % for muon conversion events with the 
momentum above a threshold momentum of 103.6 MeV/c. 

Comparing these results with the results of the reconstruction in the presence of 
the background we get the difference in tracker resolution 1.5 % and the difference 
2.7 % in overall acceptance. Therefore the tracker resolution and overall acceptance 
are not affected significantly at the considered background level. 
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